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Abstract 

We present a new scheme for extracting approximate values in "the improved perturbation method", 

which is a sort of resummation technique capable of evaluating a series outside the radius of convergence. 

Q ' We employ the distribution profile of the series that is weighted by «th-order derivatives with respect to 

o' 

^ ' the artificially introduced parameters. By those weightings the distribution becomes more sensitive to the 

i^ , "plateau" structure in which the consistency condition of the method is satisfied. The scheme works eff"ec- 
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tively even in such cases that the system involves many parameters. We also propose that this scheme has 
to be applied to each observable separately and be analyzed comprehensively. 



(N 

^ . We apply this scheme to the analysis of the IIB matrix model by the improved perturbation method 

obtained up to eighth order of perturbation in the former works. We consider here the possibility of spon- 
^ ' taneous breakdown of Lorentz symmetry, and evaluate the free energy and the anisotropy of space-time ex- 



tent. In the present analysis, we find an SO(10)-symmetric vacuum besides the S0(4)- and S0(7)-symmetric 
vacua that have been observed. It is also found that there are two distinct S0(4)-symmetric vacua that have 



jil^' almost the same value of free energy but the extent of space-time is different. From the approximate values 

1-Ch ' of free energy, we conclude that the S0(4)-symmetric vacua are most prefen-ed among those three types of 

•rH , vacua. 

>< 

. ^_ \ PACS numbers: 02.30.Mv, ll.25.-w, 11.25.Yb, 11.30.Cp, 11.30.Qc 



I. INTRODUCTION 



String theory is the unique theory that contains massless spin-two particles, i.e. gravitons UJ], 
and thus it is considered to provide an unified microscopic description of the universe including 
gravitational interactions. For this reason the string theory has been subjected to intensive studies. 
However, it is recognized that the perturbative string theory fails to single out our universe as the 
unique vacuum of the theory [j20. Therefore we are forced to pursue non-perturbative formula- 
tions. The IIB matrix model (also called the IKKT matrix model) is proposed as a constructive 
formulation of the superstring theory M,\4\- 

A significant feature of the IIB matrix model is that the space-time itself is expressed by the 
eigenvalue distribution of 10 bosonic matrices, and thus it is treated as a dynamical variable of 
the model. The origin of our four-dimensional space-time can be argued in the context of the 
IIB matrix model as a spontaneous breakdown of Lorentz symmetry. In this regard, we have to 
understand non-perturbative properties of the model. 

The mechanism of the spontaneous breakdown of Lorentz symmetry in reduced matrix mod- 
els has been examined in various approaches. It has been recognized from those works that the 
fermionic part of the action plays a crucial role O, |6J, [Zi Isi |9i llOl llin. 

Unveiling dynamical aspects of the model is, in general, quite a difficult problem. The Monte- 
Carlo method is a powerful tool for exploring such non-perturbative properties of a model. How- 
ever, it is not applicable (or at least quite difficult to apply) to the IIB matrix model due to complex 
phase of the action derived from the fermionic part' . The improved perturbation method (also 
called the Gaussian expansion method) is an alternative approach. It is considered as a sort of 
variational method dlJ, |l4i [ij, uM, livi ll8|]. It has been successfully applied to various models 
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2211 . and applications to matrix models were done in Refs. [23 



24 



m 



The application to the IIB matrix model was first achieved in Ref. ll26ll in which various patterns 
of symmetry breaking that preserve SO{d) subgroup of the original ten-dimensional rotational 
symmetry (ansatz) were examined, and as a conclusion, four-dimensional universe is the most 
preferred among them based on the comparison of free energy. The Gaussian expansion method 
was reformulated as an improvement of perturbative series expansion in Ref. [|27ll . The improved 
Taylor expansion (ITE), as it is referred, opened a way toward more general applications that 



' A novel technique called the factorization method is proposed to resolve the complex action problem in the Monte 
Carlo simulations lUuH . 



incorporate quadratic and other types of interactions. The ITE prescription was employed for the 
IIB matrix model up to fifth order of perturbation in Ref. ll27n . It has been proceeded to even higher 



orders and extended ansatz [128 . 
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30h . The mechanism of the spontaneous symmetry breakdown 



is further examined in a simplified model via the Gaussian expansion method in Ref. 13 IQ . 

In the present paper, we investigate the non-perturbative solutions of the IIB matrix model by 
this technique (which we call "the improved perturbation method" here). We focus on the possi- 
bility that the original SO(IO) symmetry of the IIB matrix model may be spontaneously broken 
to result in our universe which spreads in four directions and has S0(4) rotational symmetry^. In 
such cases it is important to see that the method is applicable to the models that exhibit phase 
transitions. In the former work ll22ll . we applied this method to the Ising model and found that the 
improved perturbation method extracts the information of the ordered phase from an expansion 
about the vacuum in the disordered phase. It is also observed that an unstable vacuum is iden- 
tified even though it has larger value of free energy than that of the stable vacuum. We expect 
that the improved perturbation method reveals diff'erent patterns of symmetry breaking that may 
be developed as unstable or metastable vacua. 

The improved perturbation method is considered as a sort of resummation of perturbative series 
by introducing the artificial parameters into the model. The approximate value of the series is 
obtained by evaluating it in the region of the artificial parameter space that realizes the principle 
of minimal sensitivity 1II6I1 . We call such a region as "plateau", in which the dependence on 
those parameters would vanish efi'ectively and the exact value would be reproduced. It works as a 
consistency condition for the parameters. 

The concept of plateau is rather obscure, and there is not yet any proper treatment of mathe- 
matical rigor. We need a practical scheme for extracting the approximate value of the series. In 
addition, it should be free from the ambiguity due to the subjectivity of recognition as much as 
possible. In the former works, mainly two approaches have been propounded. One is to take the 
values at the extrema of the improved series with respect to the artificial parameters as approximate 
values t2M- The other is to take the values which corresponds to the mode of distribution of the 



evaluated values around the accumulation of extrema. The latter is called the histogram technique 
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In this article we propose a new prescription which enables us to extract a good approxima- 



We perform the Wick rotation to the IIB matrix model and discuss in Euclidean space-time. 



tion from the improved series. It is achieved by incorporating the appropriate weightings in the 
histograms that are given by the inverse of derivatives. This method will be applicable to various 
models which have many parameters. 

It has also been customary to evaluate the various observables based on the information of 
plateau of the free energy. On extracting the approximate value of any observable, the information 
of minimal sensitivity of only improved free energy is taken into account. To be more precise, 
one takes the value of improved perturbative series of an observable at extrema of improved free 
energy with respect to the artificial parameters as approximate values of that observable. We insist 
that the improved perturbation method has to be applied independently for each observable of 
interest. 

We apply the new scheme to the analysis of the IIB matrix model by means of the improved 
perturbation method. It is based on the improved series that is obtained in former works up to 
eighth order of perturbation for SO(J)-preserving configurations as ansatz (d = 4 and 7) ll29ll . 

This paper is organized as follows. In Section UIl we provide a description of the resumma- 
tion technique called the improved perturbation method, including a review of the method. In 
Section Hill we propose a new prescription which extracts the approximate value from improved 
perturbative series. In Section|IVl we apply the improved perturbation to the IIB matrix model and 
survey the non-perturbative solutions with various patterns of symmetry breaking. Section |V] is 
devoted to the conclusion. In Appendix|Al we present a detail of analysis performed in SectionHVl 

II. IMPROVED PERTURBATION 

It is generally believed that the perturbative expansion (also including ^-expansion and s- 
expansion) is asymptotic, and diverges beyond some finite orders. Nevertheless, what we can 
evaluate for most of the theories are only those series expansions, from which we have to draw 
physical information of the theory. Therefore, we need a method to estimate the exact value from 
the series with the parameter (coupling constants, and so on) outside the original convergence 
radius. For this purpose we introduce here a method which has been successfully applied to the 
IIB matrix model and other general models. 



A. Prescription 

We implicitly use A as coupling constant and m as parameters of the model collectively such as 
masses. Let us assume that the observables of a theory would be exactly described by a function 
F(A, m). Perturbation theory provides an expansion of F as a power series of A about A = 0, with 
nth coefficient denoted as fnim): 

CO 

F(A,m) = J]rUm). (!) 

In the actual cases, we only have finite portion of the series up to order N, 

N 

F''(A,m) = J^rUm). (2) 

11=0 

The question is whether we can presume the exact value F(A, m) at the given parameters A and m 
from the series F^ above. In many cases the convergence radius for A is zero, so we can not expect 
that F^ gives a reliable approximation. 

Now we consider a modification of the series along the following prescriptions. First we per- 
form a shift of parameters: 

A^g'A, 

(3) 

m — > niQ + g''{m — itiq), 

where we have introduced ^ as a formal expansion parameter, and mo as a set of artificial param- 
eters, p and q are taken arbitrarily. We deform the series by the substitution ([3]), and then we 
reorganize the series in terms of g, drop the 0{g^^^) terms, and finally set g to 1. We obtain the 
improved perturbative series F^ as 



F^(i, m) -^ F\A, m; m,,) = F^'ign, m^ + g\m - mo)) 



(4) 



Here, we adopt a notation | ^^ ^i '^o represent the operation that we disregard the C)(g^^') terms 
and then put g to 1 . 

It should be noted that by simply setting g = 1, the modification itself becomes trivial and the 
series would be independent of the parameters mo. However, due to dropping the 0{g^^^) terms, 
the deformed series does depend on mo. 

To turn the argument around, we adopt here the principle of minimal sensitivity that the exact 
value F will be reproduced when the improved series depends least on the artificial parameters mo. 



It provides a sort of consistency condition on mo; by tuning the parameters to the solution of the 
condition, we will have a good approximation of F. 

In the above we have introduced the arbitrary parameters p, q, M. If the region emerges in the 
parameter space of mo that realizes minimal sensitivity, the approximate value must be independent 
of the choice of parameters p, q, M. However, with the limited orders of perturbation, the signal 
of minimal sensitivity is often weak in actual cases. So we have to find the optimal values of 
parameters p, q, M in order to make the signal clear. We have to keep in mind that large values of 
p and q turn to throw away the information of higher order terms of perturbative series, and taking 
larger value of M than that of N is meaningless. Therefore, we set the parameters as p = q = 1 
and M = N throughout this paper. 

B. Example 

For example, we apply the improved perturbation method to a simple 1 -dimensional function, 

FiA,m) = -^, (5) 

1 + Am 

and try to estimate F{A, m) at A = \ and m = 3/2. 

The Taylor expansion about /I = up to order A'^ gives the following series: 

N 

F^(A, m) = J^(-l T A" m" . (6) 

11=0 



This series at A'^ ^ oo has finite convergence radius \Am\ < 1 . The situation is depicted in Fig. 1 1(a) 
where F^ of various A'^ are shown as a function of m when A = I. 

Now we deform the series along the prescription in the previous section, to obtain the improved 
series F^: 

N 

F\A, m; mo) = V (-1)" {Agf (mo + g{m - mo))" „ . (7) 



n=0 



g-^.g-^l 



Here, | ^ ^ denotes disregard for the 0(g) terms followed by setting g = I. 

The improved series behaves as shown in Fig. |l(b)[ when the artificial parameter mo is taken to 
be 1.1. In this case, the functions F^{A = l,m;mo = 1.1) become close to the exact value in the 
region m ~ 1 .5, and we will have a good approximation of the original function F at m = 3/2. 

To find the optimum of parameter mo, we consider F^{A, m; mo) as a function of mo with m 
and other original parameters fixed to the specified values, (e.g. A = I and m = 3/2 in this case). 
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FIG. 1: |(a)| Taylor series of F{A, m) - j^j^ and [(b)] improved series at mo = 1.1. 



0.2 0.4 0.6 0.£ 

mo 



exact 
order 10 
order 15 
order 20 
order 25 



FIG. 2: Improved series as a function of artificial parameter rtiQ. 



It can be seen that in some region of parameter mq = 0.7 ~ 1.1, F^ stays stable and gives good 
approximation to the exact value, F = 0.4 (Fig. O. We call such a stable region as "plateau", 
where the principle of minimal sensitivity is realized and the exact value will be reproduced. 



C. Properties of improved perturbative series 



1. Details 



We will elucidate the concrete prescription of the improved perturbation method. 
The original perturbative series F^{A, m) up to order N is given in the following expression, 
where /„(m) is the nth order coefficient: 



F\A,m) = Y,X'fn{m). (8) 



n=0 



By the shift of the parameters, 



(9) 

m — > mo + g{rn — mo) , 



F^ is deformed as, 

N 

F\A, m)^Yj (^'^)" ^'('"o + g{m - mo)) 



"=" (10) 



= Z^' Z ^"(^r^(---o)^-"/r"^(mo), 

Here, each coefficient f„ is expanded in Taylor series about m = mo, with fn ^ as fcth derivative of 
/„(m) with respect to m. The series is then reorganized in powers of g. Then we drop the C)(^^^') 
terms, and finally set g to 1. As a result, we will have the improved series F^: 

N N~n , 

F'' = J]A"J]-(m- mof fP(mo) . (11) 

n=0 k=0 

It turns out that the improved perturbation method replaces each coefficient of the original 
series, fnim), by its Taylor expansion about a shifted point mo, 

N~n , 

fn(m) -^ f„im; mo) = ^ p" *^'" ~ '"^^^ frf^mo) ■ (12) 

k=o ^• 

It should also be noted that the expansion in eq. (fT2l) is taken only up to (A'^ - n)th order. This 
feature works to suppress the large ffuctuations coming from higher order coefficients, which have 
caused divergent behavior of the original series. Thus we could expect better estimates by this 
improved perturbation method. 

The important point is that the reorganized series should be truncated at finite order A'^. If we 
kept all the terms and if we could switch the order of two summations in eq. (flOl) . the procedure 
should become trivial, only to obtain the original series. 

8 



2. Convergence 

The convergence property of the improved series with respect to A is, however, not affected 
significantly. The ratio of {n + l)th coefficient against nth one, 

L^ ^ Zf=o"-' 1 (m - mo)V,S(mo) 
/„ Zf="o"n(^-^o)ViVo) ' 

has no particular structures such as singularities, at least in obvious manner. There, the parameters 

ifiQ are determined according to the following argument; minimal sensitivity condition will be 

realized ideally when pth or lower derivatives of F^ with respect to the parameters mo become 

zero: 

(jp _ 

= -FV,m;mo) 

(14) 



N p-\ 

,1=0 q=0 



f^ll , ,,. 1 



(N - n- qy. 



q ) 
We will examine the condition further in later sections. 

3. Limitations 

As mentioned before, the improved perturbation method does not much alter the convergence 
properties. There are some class of functions such that the improved series actually is convergent 
in the "plateau" region, and the procedure of improved perturbation works as a sort of analytical 
continuation. For example, the function y-^ discussed in the previous subsection belongs to this 
class rt34|l . However, most of the series that appear in physics do not have the property like this. 

If the original series is asymptotic or is evaluated at outside of convergence region, the improved 
series may also show the divergent behavior at excessively high orders. We will illustrate this 
feature by a simple example. Consider an ordinary function F, 

F{X, m) = -i— + , , ^^,, .3 , (15) 

1 + Am 1 + Am + {Amy 

and generate a Taylor series about /I = up to order A''. The radius of convergence of the original 
series is {Am) ~ 0.7. Then, we apply the improved perturbation method to this series. As shown 
in Fig. |3(a)[ the improved series dX. A = l,m = 3/2 fluctuates violently with respect to /no- The 



deviation from the exact value grows larger as the order increases. (The artificial parameter mo is 
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FIG. 3: [(a)] Large fluctuation of F as a function of mo. [(b)] Deviation between F^ and the exact value, (mo 
is chosen to be 0.95.) 

chosen to be /no = 0.95.) The improved perturbation method fails to reproduce the exact result in 
this case at high enough orders (Fig. [3(b)] ). 

In order to obtain a reliable estimate by the improved perturbation method, we have to combine 
various information of the series of different orders, and examine synthetically the whole profile 
of the series. 

4. Ring properties 

In this section, we examine how the improved perturbation method affects the function ring C°°. 
We are especially interested in the derivative operation, for the expectation value of an operator is 
related to the derivative of the free energy with the corresponding source term. 

Usually we can introduce an equivalence relation into a set of elements of function ring, and 
classify them into equivalence classes. Now we insist that the plateau criteria works for such 
classification; two elements of the function ring are equivalent if there exists intersection between 
"plateau" of each function where minimal sensitivity is realized. 

First we consider the add operation. Assume that we have two functions, F and G, and the 
improved series, F^ and G^, respectively. The procedure to obtain the improved series is basi- 
cally Taylor expansion, and thus linear. If F^ and G^ have intersecting region where minimal 
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sensitivity condition for each series is realized simultaneously, then the improved series of sum of 
the functions, (F+G) also stays stable and reproduces the exact value there. 

On the other hand, when F^ and G^ do not have overlapped parameter region, the sum (F+G) 
does not in general bear stable region in their plateaux. 

Thus, the improved perturbation method preserves the add operation among the elements of 
function ring that belong to the same equivalence class based on the presence of intersection of 
"plateau". 

Unlike the case of the add operation, there is no reliable discussion on the product operation 
of function ring. However, from several examples, it seems that the procedure also preserves the 
product operation among the elements of the same equivalence class. 

Next we consider the derivative operation. Similar to the ordinary function, the derivative of 

the improved series, F^{A, m; mo), with respect to m will be defined as follows: 

d -^Nri ^ r F^(A,m + €; ifio) - F^{A,m;mo) ..,. 

——F (A, m; mo) = lim . (16) 

dm f^o e 

We will designate as R(A, m) the region in the parameter space mo where the minimal sensitivity 
condition is satisfied. We also call it as "region of minimal sensitivity" below. 

If two regions of minimal sensitivity for physical parameters at m + e and m have overlap, 

R(i,m + 6) n R(i,m) ^0, (17) 

we would conclude from the argument for the add operation, that the derivative function defined in 
eq. (fT6l) bears the region of minimal sensitivity in the intersection of the two regions. In short, for a 
function and its derivative, the region of minimal sensitivity will be expected to emerge somehow 
close to each other. 

Otherwise, it may happen that the region of minimal sensitivity of F'^ and its derivatives are 
uncorrected. Then, in what situation does the overlap vanish, i.e. R(/l, m + e) n R(/l, m) = ? 
This question is deeply related to the critical behavior of the underlying physical model. In usual 
cases, free energy and various expectation values are regular with respect to moduli parameters, 
and the region R(/l, m) evolves regularly as well. However, at singular points of moduli space such 
as those where the symmetry of the model is enhanced, R(/l, m) should emerge in a quite diff"erent 
manner; thus even for infinitesimally small deviation of moduli parameter, R(/l, m) and R(/l, m + e) 
may not have overlaps. 

Such a situation often occurs in the physical models of our interests. For example, in Ising 
model, the region of minimal sensitivity drastically changes between high temperature phase and 

11 



low temperature phase. In the IIB matrix model (which will be examined extensively in later 
sections), we are interested in the massless region, where the parameters become singular. In these 
cases, the free energy and the expectation values of various observables show different behavior, 
and therefore we have to treat them separately. 

III. PLATEAU CONDITIONS 

In this section, we will elucidate the criteria how to determine desirable regions of artificial 
parameters, and to estimate the approximate value from the improved series. A guiding principle, 
minimal sensitivity, will be realized in regions in the parameter space where the improved function 
stays rather stable. We call such region as plateau. 

The subject of this section is to develop a concrete and objective scheme for identifying plateau. 
Besides, the notion, "stable" is obscure, and so the procedure may also provide an inductive defi- 
nition oi plateau itself. 

A. Minimal sensitivity and plateau 

The improved perturbation method introduces artificial parameters, mo, which we need to deter- 
mine by some means. Here we adopt as a guiding principle, minimal sensitivity, that the improved 
function should depend least on those parameters; this is because mo were originally introduced 
as the nominal shift of parameters, though the dependence on them appeared due to dropping the 
0{g^^^) terms. If there exists a region in the parameter space mo where the dependence on mo van- 
ishes effectively, the exact value should be reproduced there. In the following we denote simply 
by F^ the improved series F'^{A = X^,m = m*;mo) where we substitute each desirable values /I* 
and m* for A and m. This series is function of artificial parameters mo. 

The principle of minimal sensitivity is realized through the emergence of the region in the pa- 
rameter space where the improved series stays stable against any variation of artificial parameters; 
we call such a region as "plateau". What we have to do for determining mo is to identify where the 
improved series becomes stable. 

When the number of parameters are one (or at most two), the stable region of the function will 
be pointed out by just drawing a graph of it. However, the visualization will not be possible for 
the cases with more than two parameters. Moreover, even though it is possible, it may sometimes 
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lead to misjudgment by changing the scale arbitrarily. Therefore we need a concrete scheme 
for identifying the plateau applicable to multi-parameter space, which is not affected by intuitive 
guesses. 

B. Identifying plateau 

The ideal realization of the principle of minimal sensitivity may have such a property that the 
improved series is totally independent of the artificial parameters mq in a region. It involves the 
situation that all orders of derivatives of the improved series with respect to mq are zero in the 
region. However, such an ideal "plateau" is not realized in the actual cases because we have only 
finite order of series. Typical profile of the improved series that forms plateau exhibits a flat region 
in which the series fluctuates bit by bit; it would accompany a number of local maxima and minima 
there. Thus, to turn the argument around, we consider the accumulations of extrema as indications 
of plateau (or its candidates). 

It should be noted that there is also a case when the series becomes stable without forming any 
extrema in that region; we may miss such plateau by the above speculation. We will discuss this 
type of plateau in later section. 

1. Cluster identification 

Assume that we have already found extrema of the improved series in the (multi-dimensional) 
parameter space. Next issue is to identify the accumulation among them. We denote the coordinate 
of ith extrema by ^, and the distance between /th and 7th extrema by J,y. The definition of distance 
in the parameter space will be presented later. 

The distribution p(r) of the distances dij shows characteristic behavior according to the presence 
of clusters of points. Here, p(r) is given by the number of pairs (?', j) whose distance dij falls on 
to r < djj < r + 6r. As depicted in Fig. HI when one or more clusters of points exist, p(r) has a 
large peak near r = and some small bumps at larger r. Let n the number of points in the cluster. 
Then the peak height is proportional to n^, while those of other bumps are of order n. If there are 
no such clusters and the points are scattered rather uniformly, p(r) shows continuous distribution. 
This concept can be extended straightforwardly to the case in which more than one cluster are 
formed. From those distinctive properties, we are able to identify the formation of clusters and the 
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FIG. 4: Distributed points (above) and distributions of distance between points (below) for each case: (a) 



forming the cluster, and (b) straggling uniformly. 



set of points belonging to them. 

There is no intrinsic definition of distance J,y in multi-dimensional parameter space. Therefore 
we have to choose a suitable definition. It introduces metric space and must obey the triangle 
inequality. 

A naive choice is as follows: 

(18) 



, 2 _ I ^ _ ^|2 

a-- — \Xi Xj\ 



As another choice, if we can define a reasonable metric gat from some argument such as di- 
mensional analysis, we will have. 



dij= yhgds, (19) 

•j* geodesic 

along geodesic between xt and Xj. The geodesic equation may not be solved nor have any solution; 
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we would alternatively choose a linear interpolation with weights by metric as, 



dij = 2_j ^Sih) \xitk) - xitk+i)\ , xitk) = Xi + —{Xj - Xi) , 



(20) 



though this dij may not satisfy triangle inequality. 

It is not always that two extrema which close to each other in the meaning of the above distances 
belong to one stable region. This is because it is possible that between these extrema improved 
series varies its value drastically and vicinity of two extrema encounters by chance. In order to 
avert this situation, we consider one other choice of distance. This choice comes from extending 
the concept of moduli space. By the definition of distance we intend to express the degree to which 
two of the extrema of a function F'^ resemble each other. Thus we include the deviation of F^ and 
its derivatives at these two points into the notion of distance: 



d.j =\xi 



?M 



Wk 



•krN 



d'F 



dx'^ 



'kT7N 



d'^F 



^i dx'^ 



^j 



(21) 



Moreover we can also include the transition as orders of perturbation increase: 



N 
dil =\xi-Xj\+ Y^ CO MM' ^ Wk 



'kipM 



d'F 



dx^ 



krpM' 



d'F 



dx^ 



(22) 



M,M'=0 k 

The relative weights w^ and ojmm' above should be chosen according to the models considered. 

In general, we can not introduce a particular choice of metric naturally. Then we try some 
definitions of metric to a model and investigate the improved perturbative series carefully. 



2. Weighted histogram analysis 

If the improved series F^ bears a plateau, the distribution of the improved series in a region 
enclosing the plateau should form a peak corresponding to the value ofF^, which gives an estimate 
of the function on the plateau. 

Therefore, we choose a region (for ease of operation, we usually consider rectangular one) 
which encloses accumulation of extrema, and evaluate distribution p(F) of F^ ll32ll . 



P(F) 



-L 



dxl, A = {x\F <F^ix)<F+6F} 



(23) 



In multi-dimensional parameter case, the extent of the plateau is not always large nor the shape 
isometric in the parameter space. Then we need a method to extract the estimated value of F^ on 
plateau eff"ectively even when the region to be examined is taken roughly. 
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The plateau is, if ideally realized, characterized by the feature that the first and higher deriva- 
tives of F'^ should be zero or small. To enhance the contribution from such flat region, we intro- 
duce weight function to the distribution as. 



L 



p(F)= I dxw. (24) 



There may be some choices of the weight function w; for example. 



wi = ^—. , (25) 




W2 = — ^ , (26) 

dxjdxj 

and so forth. Here we denote the ?th component of coordinates by Xj. 

The histogram of weighted distribution shows sharp peak corresponding to the zeros of ^th 
derivatives of the improved series. Even though the region in question is taken roughly, the contri- 
bution from the fluctuating part of F^ is well suppressed, and we will have an estimate of F^ on 
the plateau. 

There is also another choice of weight function based on the convergent behavior of the series 
on the plateau. If the improved series converges to some value on the plateau, it implies that the 
difference between the improved series of order N and order A^ + 1 diminishes. To reflect such 
speculation, we would better choose a weight function by 

1 



W-i 



2 
pN+\ _ pN 



(27) 



In general, we combine the above weight functions on a case-by-case basis. In the following 
investigation into the IIB matrix model, we use the two weight functions vvi (|25l) and W2 (|26l) . 

3. Signal and noise 

When the number of parameters are one (or at most two), we can easily recognize the flat region 
from the graph. However, there is no such simple scheme in multi-parameter case in general. We 
are trying to identify plateau from the accumulation of the extrema as a clue. The information we 
have so far for the improved series F^ is just the location of extrema and profile of (weighted) 
distribution of F^ in some specified regions. The accumulations of extrema provides candidates 
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FIG. 5: A characteristic profile of overshoot and plateau. 



of the plateau. Among these candidates we have to distinguish the true plateaux, though yet we do 
not have definite argument on this subject. We exemplify one typical false signal, "overshoot". 

It often occurs when a discrete function is approximated by a series of polynomial (consider, 
for example, Fourier transform of rectangular waves) that the series deviates relevantly at dis- 
continuities, which remains with even higher orders taken into account. Such behavior is called 
"overshoot". In the current case, the improved series turns to be an approximation of a constant 
function by a set of polynomials because the exact value is totally independent of artificial param- 
eters. Therefore the improved series are apt to show typical overshoot behavior as well (Fig. [5]). 

The overshoot is rather isolated from other extrema that belong to a plateau; this should be 
reflected in the characteristic profile of weighted distributions. The overshoot gives minimum (or 
maximum) in its neighbor, so the histogram of the distribution weighted by vvi shows sharp peak. 



and it is cut away on lower (or upper) side (Fig. |6(a)[ ). Moreover it often happens that the peaks 
in the distributions with second and higher derivatives deviate from peaks in that with the first 
derivative. It is because the overshoot appears as isolated extrema away from those forming flat 
region. On the contrary, in the case of plateau, the histogram of the distribution weighted by wi and 
W2 shows symmetric shape and a peak of the distribution weighted by W2 lies between (or ideally 
on) peaks of the distribution weighted by wi (Fig. |6(b)[ ). For the reason stated above, a cluster of 
extrema which corresponds to overshoot can be distinguished from plateau by investigating the 
shape of distributions. 

Besides, it seems that the overshoots of different orders tend to form a line in the parameter 
space. If there is such sequential structure, we had better suspect that they are overshoots. 
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FIG. 6: (a): Distribution weighted by w\,W2 in overshoot region and (b)distribution weighted by w\,W2 in 
plateau region. 






FIG. 7: A conceptual profile of smooth plateau. 
C. Pitfalls and special cases 

I. Smooth plateau 

In the prescription thus far explained, we assumed that the improved series F^ shows fluctu- 
ating behavior on the plateau, accompanied by the accumulation of extrema. If there is a region 
where F^ becomes stable but varies gently without forming extrema (Fig. |7]), it should also be 
considered as plateau, though the above procedure based on the existence of extrema will not be 
applicable to such cases. We call this type of plateau as "smooth plateau". However, the histogram 
will show the peaked structure corresponding to this smooth plateau. The weighted distribution 
should provide a reasonable estimate of the value. 
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The problem is that we do not have any guide to locate the position of plateau in this case. One 
clue is the consistency with the plateau of other physical quantities. Though the improved series 
of different observables may show distinctive behavior, the plateau corresponding to the physical 
states should somehow be realized in each series. We have to guess the region of smooth plateau 
in parameter space based on those information. If we find a plateau in analysis of an observable, 
we must suspect the appearance of smooth plateau in neighborhood of that region in the case of 
the other observables even for absence of any extremum. 

It should be noted that we have to take care not to be confused with mere asymptotic behavior 
of the series. A concrete example will be shown in the analysis of the IIB matrix model in later 
sections. 

2. Plateau on special hypersurface 

Consider a hypersurface in the parameter space where one of the parameters is zero. If we 
change the coordinate from mo to xq = — , the first derivative with respect to the new coordinate 
becomes. 



dF^ . dF 



■N 



dxo " dmo 

Even though ^ does not vanish, ^ becomes zero at mo = 0. Thus we have to investigate the 
case separately when some of parameters take zero or oo. This is because we do not know what 
choice of the artificial parameters is preferred. 

Those cases are significant in that mo = or oo occasionally corresponds to the point where 
symmetry of the original physical model enhances; therefore the improved series may have singu- 
lar structure. 

3. Multiple plateaux 

In the multi-parameter case, it is likely that there exists more than one plateau in the parameter 
space. If each of them can be interpreted as a physical state, some correspond to stable vacua. 



while others to the unstable ones. This situation is realized in Ising model [|22|] . When we carry 
out the improved perturbation for the IIB matrix model, we expect that multiple plateaux are 
realized corresponding respectively to different vacua with different symmetry. For this reason we 
must pick up all candidates of plateau and examine these carefully. 
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D. Short summary 

In this section, we discussed the concept of "plateau" and the prescription to identify them 
in the improved series. This improved perturbation method must be applied independently for 
individual observables. 

To briefly summarize the procedure for each observable, we first find extrema of the series as 
a function of artificial parameters mo, and then identify the accumulations of the extrema which 
should be considered as candidates of plateau. Next we evaluate weighted distributions in the 
region which encloses each of the candidates, and distinguish the plateau among them. The his- 
tograms also provide estimates of the exact value from the improved series. 

This procedure assumes typical profile of the improved function. It will not be applicable when 
the function stays stable without forming extrema. In these cases, we have to locate the flat region 
referring to the information of other physical quantities, such as symmetry considerations. We 
also have to take care of the special hypersurface where some of parameters take zero or infinity; 
such cases may be relevant in the original physical models. 

IV. IIB MATRIX MODEL 

In this section we apply the improved perturbation method to the IIB matrix model and investi- 
gate its non-perturbative solutions. In particular we concentrate on the free energy and the second 
moment of eigenvalue distribution. By comparing the values of free energy of the solutions, we 
can see which solution is most preferred. In the IIB matrix model, the distribution of eigenvalues 
is interpreted as the space-time itself. Therefore we can recognize the shape of the universe by 
computing the eigenvalue distribution. To be more precise, the square root of the second moment 
of eigenvalues distribution gives the scale of extent to each direction of universe. The deviation of 
the ratio of space-time extent between directions from unity indicates the degree of anisotropy of 
the universe. 

A. Model setup 

The action of the IIB matrix model is 

5oD-YM = Tr (-^ [X,„ Xtf - ^^ [X,„ r^]] , (29) 
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where X„ (a = 1 ... 10) and ^^ (a = 1 ... 16) are all A^ x A^ hermitian matrices which belong 
to a vector and a Majorana-Weyl spinor representation of SO(IO), respectively. This action can 
be obtained by the dimensional reduction from ten-dimensional U(AO supersymmetric Yang-Mills 
action to zero dimension. It has the symmetry of all global SO(IO), \J(N) and type JIB supersym- 
metry. 

We perform the following change of scales: 

»P^iV^¥, (30) 



gl^N^A. 



Then the action takes the form: 



S =NTri-^A[X,,X,]^-^^fA^'[XaX"i']\ . (31) 

In order to carry out the perturbative expansion, we add the following propagator terms to the 
action 



5o = A^Tr 



Z', 10 1 ^° "^ 

V a,b=\ a,b,c=l ) 



(32) 



They are the most general quadratic terms. The fermionic bi-spinors are expanded by antisym- 
metric tensors by using gamma matrices. In the case of ten-dimensional Majorana-Weyl spinor, 
by considering a chirality, only gamma matrices of rank three are allowed. 

The scaling behavior of the Feynman rules for the amplitude with respect to A^ and X is given 
as follows. 

propagator for bosonic fields Xa. ~ —, 

N 

propagator for fermionic fields *F: ~ — , 

3 -point vertex: ~ N\A, 

4-point vertex: ~ NA. 

If we denote the genus of dual surface by g and the number of faces of dual surface by A^^ , we find 

that the contribution of a single diagram is 

amplitude ~ A^^d-s)^- _ (33) 

If we take the limit N ^> oo with A fixed, the planar diagrams (g = 0) give the leading contributions. 
In this case the amplitude is given by 

amplitude ~ A^^ ^^ A"f" . (34) 
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Free energy F is evaluated by the sum of connected bubble diagrams in terms of propagators m^ 
and niabc as 

The expectation value of the second momentum of eigenvalue distributions is obtained by the 
derivatives of free energy with respect to the parameters. 



dnia \ \m) 

The original IIB matrix model does not have the propagator term eq. (I32l) . Therefore we eval- 
uate the observables at m„ = niabc = 0. They are singular in this region of parameters in the above 
perturbation. We use the improved perturbation to extrapolate to the massless case. 

In general, all parameters ntat, niabc have to be taken into account. We may perform SO(IO) ro- 
tation to reduce m,,^ to the diagonal form at the stage of action. We will impose further restrictions 
from the practical reason in the later sections. 

Before going into a concrete computation, we make some assumptions. In computing the 
free energy we have to evaluate all connected bubble diagrams. As we proceed to the higher 
order of perturbation, the number of diagrams becomes enormous. We consider such a situation 
that the U(A'^) symmetry stays intact. Throughout this article we concentrate on the spontaneous 
breakdown of Lorentz symmetry^. Then we can omit the tadpole amplitude 

((X%) = (37) 

for the matrix corresponding to the adjoint representation of SUCA'^) subalgebra. By the argument 
of the gauge symmetry alone we can not forbid the tadpole amplitude of U(l) part. If it takes 
non-zero expectation value x", we perform the field redefinition 

X^^X^ + ^l^xiv. (38) 

By this redefinition we can set the tadpole amplitude to zero, 

<Tr(X")) = . (39) 



The breakdown of U(A^) gauge symmetry will be discussed in the subsequent paper 135 1 
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As we show later the effective action bears dynamical mass term. Then the zero point of free 
energy is shifted due to the redefinition of fields by eq. (1381) . Thus we can not use the zero of free 
energy as a characteristics of supersymmetry. 

Under the assumption above, the sum of connected bubble diagrams reduces to the sum of one- 
particle irreducible (IPI) bubble diagrams. Thus far we have to evaluate IPI bubble diagrams in 
order to obtain free energy. 

In the estimation of IPI bubble diagrams we can drastically simplify the computation by using 



2PI free energy [|27h . At first we compute 2PI graph amplitude Gjpi by using full propagators 



{(Xa)ij(Xb)kl) - —CabSil5jk , 

1 1 (40) 

in place of perturbative propagators lUa, niabc- Here we use double line notation, and C is the charge 
conjugation matrix. 



From G2P1 we obtain IPI bubble amplitude Fipi by Legendre transformation Bol 

10 , 10 



g(mfl, Mabc; Ca, Uabc) = G2Pl{ma, Iflabc) + J^ ^^''^" ~ 7 zli '^"*'^"«*'^ ' ^^^^ 

55 



a=\ a,b,c=\ 







dnir, 



(42) 



Ca=c(m) 



Fm = mc,=cim) ■ (43) 

Eq. (H^ corresponds to the Schwinger-Dyson equation which relates full propagators c,, and Uabc 
to perturbative propagators m^ and niabc, and eq. (|43l) produce the IPI amplitude using bare prop- 
agators m„, niabc ■ 

In the present article, we use the result of the perturbative expansion up to eighth order in A of 
Ref. \2m. 

6. Ansatz 

In this subsection, we make assumptions for the pattern of symmetry breakdown. In the case 
of the IIB matrix model, the total number of artificial parameters is quite large, namely, 10 real 
numbers for lUab (assumed to be diagonalized by SO(IO) rotation), and 120 for niabc- It will demand 



For details see also Ref. 1271]. 
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an enormous effort to search for plateau in this vast space of parameters. Therefore, we impose 
restrictions on the configuration by considering symmetry to diminish the number. In concrete 
manner, we make same assumptions for the full propagators Cab and Uabc because we compute 2PI 
bubble diagrams by using these parameters. 



In the former works 1126 . 



m 
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2911 . the configurations called SOid) ansatz have been inten- 
sively examined which preserve SO(J) rotational symmetry. The guideline of choice is described 
as follows. First, SO(J) subgroup of SO(IO) is chosen to which directions the expectation values of 
fermionic two-point function Uahc are zero, d is taken from 1 to 9. Toward the rest of the directions 
Uabc may have non-zero value. Since Uabc is a rank three anti-symmetric tensor, a single non-zero 
component oiuabc brings out three-dimensional subspace by permutation of indices. Thus, (10- J) 
dimensional part would naturally be decomposed into multiples of three-dimensional blocks. Fur- 
thermore, those blocks are subjected to the permutation symmetry under the interchange with each 
other. In this way, SO(J) {d = 1 • • • 7) symmetric ansatz has been investigated in former works. 

Among those choices shown above, we examine d = A and d = 1 cases in particular. It is 
reported in Ref. [27] that d = 5 and 6 cases reduce to S0(7) ansatz, while d = 2 and 3 cases to 
S0(4) ansatz. J = 1 had no solution. 

The preserved symmetry and the explicit forms of the full propagators for d = A and d = 1 
cases are given as follows. 
SO(7) ansatz: S0(7) x S0(3) 

Cab = diag( 7crs, 3c2's ) , ^ = mF'-^''^ , (44) 

SO(4) ansatz: S0(4) x S0(3) x S0(3) x Z2 

Cab = diag( 4crs, 6c2's ) , ^ = ^(F^-^-^ + F^-'^'i") , (45) 

V2 

the Z2 factor stands for the permutation symmetry between two S0(3) factors. 

For the ansatz presented above, we evaluate the free energy given by the sum of IPI bubble 
diagrams. It is obtained as a series of the coupling constant A whose coefficients are functions of 
the parameters mi, mj, and m. Next, we calculate the second moment of eigenvalue distribution 
by differentiating the free energy with respect to mi and m2. Finally, we apply the improved 
perturbation method to these perturbative series: 

F^{A,mi,m2,u) — > F^(i,(mo)i,(mo)2, "o) 

= F\gA, {mo} + g{{m] - {mo})) „ , (46) 
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{m} and {mo} represent collectively the set of parameters {mi , m2, m} and {(mo)i , (mo)2, mo}, respec- 
tively. We put mi = m2 = m = so that it describes the original IIB matrix model. At this stage, 
we set /i = 1 . It is because the coupling constant A can be absorbed by the redefinition of the fields. 
We introduce the new variables xi as 

Xl = - — -, X2 = - — -, ^3 = — , (47) 

(mo)i (mo)2 mo 

and use these parameters in the following section. 

C. Results in S0(4) ansatz 

In this subsection, we present the results in the case of S0(4) ansatz. We apply the improved 
perturbation method to the perturbative series for the free energy, the second moment of eigenvalue 
distribution for X" (a = 1 • • • 4), which is denoted by ci, and that for X" (a = 5 ■ ■ ■ 10), which is 
denoted by C2. 

As explained in Sec. HIl we first search for the extrema of these functions with respect to the 
artificial parameters x,. Then, we find the accumulations of extrema and consider them as the 
candidates of plateau. We compute the weighted distribution of the functions in those regions. 
We use here the first derivatives (|25l) and the second derivatives (|26l) in particular. Finally, we 
identify the region as plateau in which two distributions overlap with each other. The graphs of 
distributions in individual regions are shown in Appendix lAl 

1. Free energy 

Fig. [8] shows the distribution of extrema of the improved series of free energy. As is seen 
in Fig. m there are several accumulations of extrema. We pick up some accumulations as the 
candidates of plateaux and compute the distributions with weight functions wi and wj- Besides 
most of accumulations show the hopeless distributions which spread randomly, two regions A and 
B which are depicted in Fig. [8] have a hopeful distribution. In the following, we refer to only 
hopeful candidates. 

Two regions A and B are hopeful candidates of plateau, although the behavior of the function 
of the seventh order in region A is rather unstable. The similar situation also occurred in the 



previous work [|29h in which the positions of extrema show a different pattern for the seventh order 
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FIG. 8: Distribution of extrema of improved free energy for S0(4) ansatz plotted in X1-X2 plane. 

of perturbation. In the region A we read the value of free energy to be -1 ~ 2. And in the region 
B, we also read the value of free energy to be -1 ~ 2. 

From this accordance between region A and B about the value of free energy, we may expect 
that two regions belong to one plateau. However from the examination of the other observables in 
the later subsections, we conclude that two regions are distinguished each others and there are two 
plateaux. 

2. ci 

Fig. m shows the distribution of extrema of the improved series of ci . As mentioned in former 
section, although c\ is the derivative function of free energy with respect to mi, we have to apply 
the improved perturbation separately. It is because the behavior of the function becomes singular 
in the region of moduli space that corresponds to the massless case. It is seen by comparing Fig. [8] 
with Fig. m Two distributions of extrema show the different pattern. 

Among accumulations, the regions A, B and C are promised candidates of plateau. The distri- 
bution in the region A shows a strange behavior. As the order of perturbation increases, the value 
at the peak of distribution shifts to the right. When we approximate an infinitely large value by 
a finite series, such a behavior may often be observed. Therefore this property indicates that the 
exact value is quite large. Up to the eighth order of perturbation, we conclude that ci takes value of 
40 ~ 60. From the distribution in the region B, we read the value of ci to be 1.5 ~ 2.1. Although 
regions A and B point the similar value of free energy, these regions shows clearly different distri- 
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FIG. 9: Distribution of extrema of the improved ci for S0(4) ansatz. 
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FIG. 10: Distribution of extrema of the improved C2 for S0(4) ansatz. 

butions with respect to C]. Therefore we conclude that these regions belong to different plateaux. 
On the other hand, from the region C, we read the value of ci to be 0.25 ~ 0.26. Thus we conclude 
that there are three plateaux. In particular, it is interesting that the region C shows the same value 
as the region C of C2 and S0(7) ansatz which will be shown in later subsections. The significance 
of this accordance will be discussed later in Subsection IIVEI 

3. C2 

Fig. [TOl shows the distribution of extrema of the improved series of C2. We consider the regions 
A, B, C and D as candidates of plateau. From the distribution in the region A, we read the value of 
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FIG. 1 1 : Distribution of extrema of the improved free energy for S0(7) ansatz plotted in xi -xj plane. 

C2 to be 0.0 ~ 0.2. From the distribution in the region B, we read the value of C2 to be 0.13 ~ 0.15. 
From the distribution in the region C, we read the value of C2 to be 0.25 ~ 0.26. 

In the similar fashion, region D indicates that C2 takes the value 0.21 ~ 0.25. From the similarity 
of the distributions in the region C and D, we conclude that these two regions belong to one plateau. 
Thus we choose the region C as a representative of this plateau. 

D. Results in S0(7) ansatz 

In this subsection, we show the results for the case of S0(7) ansatz. The method of analysis is 
the same as that in the S0(4) case. The graphs of distributions in individual regions are shown in 
Appendix El 

I. Free energy 

Fig. E] shows the distribution of extrema of the improved series of free energy. Though it 
seems that the region B contains three distinct accumulations, we treat them as one region. It is 
because the distributions of three regions show the similar shape, then it is indicated that all three 
accumulations belong to one structure. Thus, there are two candidates of plateau. 

From the distribution in the region A, we read the value to be 5 ~ 7. From the distribution 
in the region B, we read the value of the free energy to be -10 ~ 0. The extrema in the region 
B are considered to be overshoots, because the weighted distribution with the second derivative 
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FIG. 12: Distribution of extrema of improved ci for S0(7) ansatz. 

deviates from that with the first derivative, and it is a characteristic feature of overshoot. Moreover, 
the behavior of the improved series of other observables around this region looks despairing. We 
conclude that it is indeed the overshoot. 

2. ci 

Fig. [12] shows the distribution of extrema of the improved series of Ci. Regions A and C 
are hopeful accumulations. From the distribution in the region A, we read the value of ci to be 
0.7 ~ 0.8. From that in the region C, we read the value of ci to be 0.25 ~ 0.26. 



3. C2 



Fig. [13] shows the distribution of extrema of the improved series of C2. Regions A and C are 
hopeful regions. From the distribution in the region A, we read the value of C2 to be 0.08 ~ 0.12. 
From that in the region C, we take the value of C2 to be 0.25 ~ 0.26. 



E. SO(10)-symmetric vacuum 

From the investigation of the observables ci and C2, we find regions in both S0(4) and S0(7) 
ansatz that have an interesting property. It is shown by region C in Figs. [9] and [10] for S0(4) ansatz, 
and in Figs. [12] and [I3] for S0(7) ansatz as well. The approximate values of Ci and C2 estimated 
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FIG. 13: Distribution of extrema of improved C2 for S0(7) ansatz. 
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FIG. 14: Dependence of the improved free energy with respect to x\ and X2 at X3 - 0.1 
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dependence of the improved free energy at xi = ^2 = 0.15 (b) for S0(4) ansatz 



from the distributions are 0.25 ~ 0.26, and they are almost the same. In addition, the region C 
is located at xt, ~ 0, and thus it implies that the bi-fermionic variable mi vanishes. Therefore, we 
identify this region as the plateau that corresponds to the SO(10)-symmetric vacuum. 

It may seem strange that there is no such accumulation of extrema found for the improved 
free energy around region C above. In fact, this region appears as "smooth plateau" discussed in 
Sec. Uni By plotting the dependence of the improved free energy (Fig. O, we can recognize that 
the improved free energy stays almost constant in region C, and thus the minimal sensitivity is 
indeed realized there. The approximate value of the free energy is obtained from the distribution 
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shown in Appendix lAl as 6 ~ 7. 

The above observations are obtained for both S0(4) ansatz and S0(7) ansatz, and the approx- 
imate values of the free energy, ci, and cj are the same between those two ansatz. It is consistent 
with the speculation that these solutions obtained for two distinct ansatz are actually identical, in 
which the SO(10)-symmetry is realized. 

F. Results in IIB matrix model 

In the examinations of three observables, free energy, c\ and C2, for the S0(4) ansatz we find 
three sets of candidate of plateau. 

First plateau (we call this region as region A) indicates that the value of free energy is -1 ~ 2, 
that of c\ reaches 50 and that of ci is less than 0.2 up to the eighth order of perturbation. This 
plateau corresponds to an anisotropic vacuum that preserves S0(4) subgroup. The ratio of extent 
of the eigenvalue distributions between the 4-dimensional part and the 6-dimensional part attains 
at 20 at the stage of the eighth order of perturbation. The behavior of distributions as increase of 
a order of perturbation indicates that the exact value of ratio is quite large. This solution has a 
property that the degree of anisotropy is large (probably infinite). Next plateau (we call this region 
as region B) indicates that the value of free energy is -1 ~ 2 and that of C\ is 1.5 ~ 2.1 and that 
of C2 is 0.13 ~ 0.15. This plateau also corresponds to an anisotropic vacuum that preserves S0(4) 
subgroup. The ratio of extent of the eigenvalue distributions between the 4-dimensional part and 
the 6-dimensional part is obtained by 3.1 ~ 4.0. This solution has a property that the degree of 
anisotropy is finite. 

The third candidate appears in a different manner. This is the plateau corresponding to a 
S0( 10)- symmetric vacuum which was discussed in the previous subsection. For the observables 
C\ and C2, it indicates that both of ci and C2 take the value 0.25 ~ 0.26. For the free energy this 
plateau appears as "smooth plateau". By computing the distribution of free energy we conclude 
that the approximate value of free energy is 6 ~ 7. The graphs of distributions in this regions are 
shown in Appendix lAl 

Similar argument can be applied to the case of S0(7) ansatz. There are two plateaux; one of 
which indicates that the value of free energy is 5 ~ 7 and that of ci is 0.7 ~ 0.8 and that of ci is 
0.08 ~ 0.12, and it corresponds to an anisotropic vacuum that preserves S0(7) subgroup. The ratio 
of the extent between the 7-dimensional part and the remaining 3-dimensional part is 2.4 ~ 3.1. 
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The other plateau corresponds to a SO(10)-symmetric vacuum. This plateau indicates that the 
value of free energy is 6 ~ 7 and that of ci and cj is 0.25 ~ 0.26. Those values are the same as 
those obtained for the third plateau of S0(4) ansatz. It is consistent with the speculation that they 
are actually identical, and SO(IO) symmetry is realized. 

V. CONCLUSION 

We have investigated the properties of the improved perturbation method in much detail. This 
method is a sort of systematic optimization of a series expansion obtained by the perturbation the- 
ory, achieved by reorganizing the series with the artificially introduced parameters. In particular, 
we focus on the realization of the principle of minimal sensitivity, which provides a condition for 
the parameters that the improved series should be least dependent on them. A region in the param- 
eter space that satisfies the above condition is called as "plateau". The exact value of the series 
would be reproduced on the plateau. Thus the efficient scheme for the identification of the plateau 
and evaluation of the approximate value becomes a major issue in the application of the improved 
perturbation method. 

We proposed a new scheme for evaluating the improved series which yields good approximate 
values. It relies on the distribution profile of the series that is weighted by the first and the second 
derivatives with respect to the artificial parameters. By incorporating those weightings, the distri- 
bution becomes more sensitive to such a behavior that the series stays almost constant in a certain 
region of the parameters, i.e. the plateau condition be satisfied. It is also noted that this scheme 
works effectively even when there are a number of parameters. 

We insist that the improved perturbation method should be applied independently for each indi- 
vidual observable of interest. It is because the prescription of improvement is quite nonlinear and 
the formation of plateau may occur in a different manner especially at low orders of perturbation. 
If they are supposed to coincide at high enough orders, the information of plateau of other observ- 
ables may help to clarify the behavior of the improved series. Thus those information should be 
referred synthetically. 

We studied the IIB matrix model by the improved perturbation method in order to examine 
the non-perturbative solutions of the model. We concentrate on the possibility of the spontaneous 
breakdown of Lorentz symmetry in the present article. We applied our scheme to the analysis 
of the improved series up to eighth order of perturbation for S0(4) and S0(7) ansatz that have 
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been obtained in former works. As a result, we found a new solution that corresponds to the ten- 
dimensional universe besides the four- and seven-dimensional solutions already found previously. 
A key to this detection resides in that we performed the analysis independently for individual 
observables and combined the information on the plateau of them. 

The existence of the plateau of the SO(10)-symmetric vacuum is quite substantial because it 
provides an evidence for the validity of application of the improved perturbation method to the IIB 
matrix model. It ensures us that the investigation encloses the whole moduli space of the model 
including the SO(IO) symmetry and that the improved series indeed describes the original SO(IO)- 
symmetric IIB matrix model. It is a non-trivial issue when the model exhibits phase transitions, as 



was discussed in the application to the Ising model 12211 . 

Based on our new scheme for the improved perturbation method, we found that the value of 
free energy becomes large in the order of four-, seven- and ten-dimensional universe. Thus we 
conclude that the four-dimensional universe is more preferred than seven- and ten-dimensional 
universe. 

It is interesting that there are two distinct solutions that correspond to four-dimensional uni- 
verse. One is the solution in which the degree of anisotropy is quite large, and the other is the 
solution in which the degree of anisotropy is finite. The estimated values of free energy are al- 
most the same for those two solutions. It may suggest an interesting possibility that the inflation 
is caused by the phase transition from one S0(4) vacuum with finite degree of anisotropy to the 
other S0(4) vacuum that has infinitely large extent in four directions, i.e. our universe. 
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APPENDIX A: DETAILS OF ANALYSIS IN IIB MATRIX MODEL 

In this appendix we show the graphs of distributions in individual ansatz and these regions 
which are candidates of plateaux. 

As we can see in Figl8lfT3l for each observables there are many accumulations of extrema. We 
must compute the distribution with weight function wi, Wi for each region. Most of regions show 
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hopeless distributions. Then we select only the hopeful regions here and depict the graphs of the 
distributions in these regions. The locations of these regions in the artificial parameter space are 
depicted in Fig[ 
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2. SO(7) ansatz 
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3. SO(10)-symmetric vacua 
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